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ABSTRACT 

We present velocity- resolved spectro-astrometric imaging of the 4.7 /im rovibrational lines of CO 
gas in protoplanetary disks using the CRIRES high resolution infrared spectrometer on the Very 
Large Telescope (VLT). The method as applied to three disks with known dust gaps or inner holes 
out to 4-45 AU (SR 21, HD 135344B and TW Hya) achieves an unprecedented spatial resolution of 
0.1 — 0.5 AU. While one possible gap formation mechanism is dynamical clearing by giant planets, 
other equally good explanations (stellar companions, grain growth, photo-evaporation) exist. One 
way of distinguishing between different scenarios is the presence and distribution of gas inside the 
dust gaps. Keplerian disk models are fit to the spectro-astrometric position- velocity curves to derive 
geometrical parameters of the molecular gas. We determine the position angles and inclinations of 
the inner disks with accuracies as good as 1-2°, as well as the radial extent of the gas emission. 
Molecular gas is detected well inside the dust gaps in all three disks. The gas emission extends to 
within a radius of 0.5 AU for HD 135344B and to 0.1 AU for TW Hya, supporting partial clearing 
by a < 1 — 10 Mjup planetary body as the cause of the observed dust gaps, or removal of the dust 
by extensive grain coagulation and planetesimal formation. The molecular gas emission in SR 21 
appears to be truncated within ~ 7AU, which may be caused by complete dynamical clearing by 
a more massive companion. We find a smaller inclination angle of the inner disk of TW Hya than 
that determined for the outer disk, suggestive of a disk warp. We also detect significant azimuthal 
asymmetries in the SR 21 and HD 135344B inner disks. 

Subject headings: techniques: high angular resolution - planetary systems: protoplanetary disks - 
stars: individual (EM* SR 21, HD 135344B, TW Hya) - ISM: molecules 
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1. INTRODUCTION 

The recent identification of a significant class of proto- 
planetary disks with gaps and inner holes in their ra- 
dial dust distribution has given rise to a debate regard- 
ing their origin (Najita et al. 2007). Dynamical clearing 
by young giant planets is the most popular explanation 
(Strom et al. 1989; Skrutskie et al. 1990), but other po- 
tential explanations exist. The presence of an inner re- 
gion in the disk that is optically thin to continuum pho- 
tons from near-infrared to millimeter wavelengths is not 
sufficient evidence for a newly formed planet, but simply 
indicates a deficit of small dust grains. 

A stellar companion orbiting at >1AU, but too close 
to the primary to be resolved with direct imaging, would 
clear out a large hole completely (Artymowicz & Lubow 
1994), such as is the case for CoKu Tau/4 (Ireland & 
Kraus 2008). However, it has been demonstrated that a 
companion body smaller than a certain threshold mass 
will allow some gas and small dust grains to accrete 
through its orbit (Lubow ct al. 1999; Kley et al. 2001), 
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specifically if the Roche lobe of the companion is smaller 
than the disk scale height. The threshold mass depends 
on many parameters, but reported values are between 1- 
10 Jupiter masses (Lubow et al. 1999), and the presence 
of gas inside a dust hole or gap therefore seems to rule 
out clearing by a stellar companion. Photo-evaporation 
will create holes in the radial gas and dust distribution 
extending all the way to the stellar surface, if the disk 
is not too massive (Alexander & Armitage 2007). Fi- 
nally, removal of dust opacity by a systemic growth of 
dust particles to meter-sized bodies, or even planetes- 
imals, will remove the dust emission, but preserve the 
gas content (DuUemond & Dominik 2005). The various 
scenarios may be distinguished by the radial distribution 
of molecular gas inside the optically thin regions, since 
some of the mechanisms for gap formation tend to re- 
move or destroy molecular gas over a specific range of 
radii, while others do not. 

In this paper we distinguish between "gaps" , which are 
radial clearings leaving an inner disk intact, and "holes", 
which refer to cases where no inner disk is present. 
We also distinguish between "dust opacity gaps/holes" 
(shortened to dust gaps) and "gas gaps/holes", each re- 
ferring to separate components traced by different sets of 
observational data. Specifically, in this paper we explore 
whether (molecular) gas gaps are present in disks with 
dust gaps. 

Regardless of the mechanism forming actual disk gaps 
and holes, it also essential to measure the radial gas sur- 
face density to understand the evolution of solids and the 
formation of planetesimals. The dynamics of large par- 
ticles, in particular of radii 1-100 cm, are intimately cou- 
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Fig. 1. — A small section of the spectral range covered for SR 
21 showing the flux spectrum (top), the atmospheric transmis- 
sion spectrum {middle) and the spectro-astrometric signal (bot- 
tom). The units on the vertical axis are in pixels, relevant for the 
spectro-astrometric signal; the other two curves are scaled for clar- 
ity. The emission lines in the j; = 1 — rovibrational band (unless 
otherwise indicated) intrinsic to the SR 21 disk are marked with 
vertical dashed lines. 

pled to the radial profile of the gas surface density. This 
is because 10 — 100 cm sized solids tend to migrate ra- 
dially against the gas pressure gradient (Weidenschilling 
1977). Thus, the presence of "pressure bumps" in the 
disk may concentrate the densities of such "pebbles" and 
"boulders" , potentially leading to the formation of plan- 
etesimals (Johansen et al. 2007). Further, gaps in the ra- 
dial gas distribution will block the inwards migration of 
particles rich in volatiles, strongly affecting the chemistry 
of the inner disk (Ciesla & Cuzzi 2006). These processes 
take place in the planet-forming zones of protoplanctary 
disks (< 10 AU) requiring imaging of (molecular) gas on 
1-10 milli-arcsec scales. 

In this paper we introduce spectro-astrometric imaging 
of molecular gas in protoplanetary disks as a technique 
that allows a spatial resolution better than 1 milli-arcsec 
(~ 0.1 AU at the distance of the nearest star forming 
regions). The data not only confirm that molecular gas 
in Keplerian motion is present well inside the dust gaps 
of three disks, but also measure the distribution (radially 
and azimuthally) of the emitting gas relative to the dust. 
The disks, TW Hya, SR 21 and HD 135344B, are part 
of the class of so-called "transitional" or "cold" disks, 
which show dust gaps or holes out to 4, 18 and 45 AU, 
respectively (Calvet et al. 2002; Brown et al. 2007), at 
distances of 51, 120 and 84 pc (Mamajek 2005; Loinard 
et al. 2008; Sylvester et al. 1996). HD 135344B and SR 21 
are examples of disks with dust gaps, since strong near- 
infrared veiling requires the presence of an inner disk 
component at < 0.5 AU (Brown et al. 2007). TW Hya 
appears to have a 4 AU dust hole (Calvet et al. 2002), al- 
though recent near- and mid-infrared interferometry has 
detected some dust emission within 4 AU (Eisner et al. 
2006; Ratzka et al. 2007). 

While HD 135344B and TW Hya both show varying 
degrees of accretion at relatively low rates (5 x 10~^'^ to 
5 X lO^'' Af0 yr~^), inferring the presence of some gas 
very close to the central stars (< 0.1 AU) (Muzcrolle 
et al. 2000; Garcia Lopez et al. 2006), this does not trace 
the distribution of gas in the disk at larger radii. No 
estimate of the accretion rate in SR 21 is available, but 
the weak Pf/3 line observed at 4.65 y^m suggests that it 
is low. The mere presence of rovibrational CO lines sug- 



gests that the gas is relatively close to the central star 
since high temperatures are required to excite them (Na- 
jita et al. 2003; Salyk et al. 2007). However, the CO 
excitation temperature itself is not enough to locate the 
gas since it is known that rovibrational lines are often 
non-thermally excited to large radial distances (Blake & 
Boogert 2004; Goto et al. 2006), introducing significant 
ambiguities in conclusions based on an assumption of 
thermalized lines. 

2. OBSERVATIONS 

Traditional infrared imaging cannot resolve the struc- 
tures in the inner planet- forming regions of circumstellar 
disks. Even with modern adaptive-optics assisted tele- 
scopes, the very best resolution achieved at 2-5 /xm is 
roughly defined by the airy ring at 1.22A/I? = 60 — 150 
milli-arcseconds for an 8.2 m telescope, corresponding to 
6-15 AU at the relevant distance of 100 pc (Goto et al. 
2006). Infrared interferometers operating at 2 /zm can 
reach resolutions about 20-30 times better, but cannot 
trace the gas kinematics at 1-10 AU in protoplanctary 
disks. 

Here, we employ spectro-astrometry as a technique 
for obtaining very high spatial and spectral resolution 
line imaging. This technique has been applied at visi- 
ble wavelengths (Takami et al. 2001; Porter ct al. 2005), 
but not yet to longer infrared wavelengths and to molec- 
ular lines in protoplanetary disks. We demonstrate the 
feasibility of using spectro-astrometry for imaging the 
Av = 1 rovibrational band of CO at 4.7 /im to a preci- 
sion <1 milli-arcsec. The method is based on the fact 
that the centroid of the image of a point source can be 
determined to much higher precision than the size of the 
point spread function (PSF), by several orders of mag- 
nitude. By measuring the centroid (1st order moment) 
in the spatial direction of a spectrally resolved emission 
line, one can construct position- velocity diagrams, given 
a certain position angle of the entrance slit of the spec- 
trometer on the sky. 

For this purpose we used the recently commissioned 
high resolution (A/AA ~ 100, 000) infrared spectrometer, 
CRIRES, mounted on UTl of the Very Large Telescope 
of the European Southern Observatory (Kaeufl et al. 
2004). The instrument is fed by the Multi- Application 
Curvature Adaptive Optics (MACAO) system, making 
it ideal for spectro-astrometry. The spectro-astrometric 
technique as applied to a larger suite of CRIRES spec- 
tra is described in greater detail in an upcoming paper 
(Pontoppidan et al., in prep.). 

We obtained spectra centered at 4.715 ^m in the CO 
P-branch at 6 different position angles (PA) for each 
disk. The slit width was 0'.'2, compared to the diffraction- 
limited PSF core of 0'.'18. The 2-D spectra were pro- 
cessed using standard reduction techniques for infrared 
spectroscopy. The formal spatial centroid of the 2-D 
spectra is calculated as: 

Xiv) = c^^^^P-^^, [pixels] (1) 

where Xi {v)—Xo is the center of pixel i relative to the con- 
tinuum centroid at velocity v in the spatial direction and 
Fi (v) is the flux in that pixel. The sum is over a given vir- 
tual aperture in the spatial direction, in our case ^ ±20- 
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TABLE 1 
Log of observations 



Star PA Obs. Date Int. Time Spectral Range 



SR 21 


190° 


30/8/2007 32 


min 


4.660-4.770 (Um 


SR 21 


250° 


30/8/2007 32 


min 


4.660-4.770 ixm 


SR 21 


310° 


31/8/2007 32 


min 


4.660-4.770 fiia 


HD 135344B 


180° 


22/4/2007 20 


min 


4.645-4.755 |Um 


HD 135344B 


240° 


4/9/2007 20 


min 


4.660-4.770 |Um 


HD 135344B 


300° 


5/9/2007 20 


min 


4.660-4.770 ^tm 


TW Hya 


243° 


26/4/2007 40 


min 


4.660-4.770 fim 


TW Hya 


333° 


26/4/2007 40 


min 


4.660-4.770 fim 



of the PSF, or 400 milli-arcsec. C is a correction factor 
of order unity that takes into account that not all of the 
source flux is enclosed in the virtual aperture. For our 
choice of aperture, which maximizes the signal-to-noise, 
C = 1.3. This value was confirmed by measuring the am- 
plitude of the spectro-astrometric signature as a function 
of virtual aperture size. Further, the total centroid is di- 
luted by the continuum flux such that the true centroid 
of the line is Xi{v) = X{v){l+F^{v)/Fi{v)), where FJFi 
is the total continuum-to-linc ratio at velocity v. SR 21 
and HD 135344B have Fc/Fi ^ 5, lowering the spectro- 
astrometric sensitivity by a factor 6, while TW Hya has 
much smaller continuum dilution with Fc/Fi ~ 0.7. 

The PAs are paired such that each spectrum has a 
counter-spectrum obtained at an angle rotated by 180° 
(Brannigan et al. 2006). By pairwisc subtracting the 
parallel and anti-parallel spectro-astrometric signals, ar- 
tifacts stemming from instrumental effects are removed 
while any real signal is retained, i.e. X{v) = (X^o — 
Xjj^qo)/2. Grotmd based 4.7/im spectra arc filled with 
numerous strong absorption lines due to atmospheric 
molecules, including CO itself. We use these strong spec- 
tral features to test for the presence of artifacts in our 
spectro-astrometry. For the data presented here, no tel- 
luric features produce false spectro-astrometric signals 
in excess of the statistical errors after pair subtraction 
(down to X{v) ^ 200 — 500 /iarcsec), ruling out the pres- 
ence of artifacts as discussed in Brannigan et al. (2006). 
This is illustrated in Figure 1, which shows that the 
spectro-astrometry is consistent for all the lines intrin- 
sic to the disk, but is not correlated with the shape of 
the telluric absorption spectrum. 

The CRIRES setting employed covers ~8 CO t;=l-0 
lines and ~13 CO ?;=2-l lines. The spectro-astrometric 
lines from each vibrational band are co-added to reveal 
position-velocity diagrams at 3 position angles for each 
disk (2 for TW Hya), obtained with relative position an- 
gle spacings of 60° (90° for TW Hya). Care was taken 
to avoid lines contaminated by other transitions or fore- 
ground absorption. In particular the three lowest J tran- 
sitions for SR 21 show an absorption component due to 
gas in the foreground molecular cloud. While such ab- 
sorption will not affect the spectro-astrometry, it will af- 
fect the line flux proflle. Thus these lines were not used 
for the co-added flux spectra and spectro-astrometry. 
The flux spectra were obtained by division with an early- 
type telluric standard (HR 6084 for HD 135344B and SR 
21; HR 4023 for TW Hya). 

Ideally, each position angle for a given source should 
be obsc^rved within a time span of a few days. This is 
the case for the observations presented, except for the 
PA = 180° of HD135344B which was observed 4 months 



prior to the other position angles of that source. The 
observations are summarized in Table 1. 

3. RESULTS 

3.1. Keplerian disk models 

An emission line originating in a Keplerian disk is ex- 
pected to show an antisymmetric position- velocity spec- 
trum, with the blue- and red-shifted sides of the line 
offset in opposite directions. The CO position-velocity 
diagrams for SR 21, HD 135344B and TW Hya are shown 
in Fig. 2. The CO A?;=l lines typically probe molecu- 
lar gas at temperatures of 500-1200 K (Blake & Boogert 
2004). It was possible to extract high quality spectro- 
astrometric signals from both CO v =1-0 and CO v=2-l 
lines in SR 21. All the astrometric CO spectra, with the 
exception of TW Hya at PA=63°, show clear Keplerian 
signatures. 

In order to infer the physical parameters of the gas 
emission, a simple model of a Keplerian disk is fitted 
to each set of position- velocity diagrams. The model is 
constructed as a discrete set of emitting rings with a ra- 
dial temperature dependence determined by a power law: 
T = Ti Au(-R/1 AU)^^, and a constant radial surface den- 
sity. The free parameters in the model are the radius of 
the innermost ring, the PA of the major axis of the pro- 
jection of the disk on the sky, the inclination of the disk 
and the stellar mass. The models were fitted simulta- 
neously to all the spectro-astrometric position angles, as 
well as the line flux spectrum. 

The spectro-astrometric signature is quite sensitive to 
the PA and the weighted average of the radius of the 
emitting gas, and these parameters are generally deter- 
mined to high accuracy. The stellar mass and disk incli- 
nation are degenerate, tracing a narrow strip of good fits 
as seen in the surfaces shown in Fig. 3. While the 
spectro-astrometry measures the radius of the weighted 
average of the line emission, the constraints of the cor- 
responding line profile add information about the radial 
extent of the gas. 

For SR 21, the data could be flt with a single ring 
{AR/R ^ 1) at 7-7.5 AU because the line shape is double 
peaked (no emission from large radii) and has steep wings 
(no emission from small radii). Attempts at fitting the 
SR 21 data with a temperature distribution lead to very 
similar results. Essentially, it is required that the line 
surface brightness drops rapidly with increasing radius 
so as to not contribute to the total line profile. However, 
an actual truncation of the gas at some outer radii is not 
strictly required, since it is also possible that it is simply 
not warm enough to be excited. This scenario constrains 
the temperature profile to q > 0.4. We consider the 
more complex model with relatively cold gas extending 
outwards from an inner cut-off radius to be more physical 
and therefore adopt it, even though a single ring also fits 
the data. 

For HD 135344B and TW Hya, a ring was insufficient 

to fit both the spectro-astrometry and the line profile in 
detail, suggesting that the line emission comes from a 
range of radii (see Table 2). The cause of this is most 
clear for HD 135344B, which has a line profile with broad 
wings. For TW Hya, we use the temperature profile es- 
timated for the dust using near-infrared interferometry, 
TiAu = 360 K and q = 0.4 (Eisner et al. 2006). This 
value for q corresponds to optically thin dust in LTE, 
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Fig. 2. — Top panel: average CO P-branch lines in the v=l-0 fundamental rovibrational band at 4.7 /xm for the three cold disks. For 
SR 21, the 11=2-1 lines are displayed as well. Bottom panels: Spectro-astrometry for the same lines as defined in eq. 1; the y-axis offset is 
for the line+continuum signal with respect to the continuum emission centroid. The red curves show the best fitting disk models. At the 
bottom of each spectrum the transitions used for the composite are indicated. The telluric CO lines are visible as a gap in the spectra. 
The errors have been propagated from the pixel RMS variations in the 2-D spectrum, under the assumption that the error is dominated 
by background noise (a good approximation at 4.7/im). 
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Fig. 3. — surfaces for the best fitting Keplerian disk models. The contours are at 68, 95 and 99% confidence levels, corresponding to 
1, 2 and 3cr. 



(7 = 2/(4 + where /? = 1 is the dust opacity power law 
index for small (non-grey) dust grains. For HD 135344B, 
Ti AU = 790 K and 9 = 0.4 provide a reasonable fit to the 
line profile, although the line wings are so broad that a 
single power law temperature profile fails to provide a 
perfect fit. A shallower temperature profile with q < 0.4 
over a range of radii may improve the fit, but explor- 
ing more complex models requires the inclusion of more 
physics, such as fluorescence excitation, in the excitation 
temperature profile, and is beyond the scope of this pa- 
per. Further, the line is asymmetric and the blue line 
wing is stronger than what can be fitted by our simple 



model. For SR 21, Tvau = 300 K and q = 0.4 was as- 
sumed, relevant for LTE conditions. Note, however, that 
given the strength of the CO v = 2 — 1 lines, the vibra- 
tional ladder is probably non-thermally excited, even if 
the rotational is not. 

The most obvious result is that the molecular emission 
originates well within the dust gaps reported in the liter- 
ature, and that the radial distribution of the line emission 
is not the same for the three disks. SR 21 appears to be 
devoid of molecular gas within ~ 5 AU, while the molec- 
ular gas extends to within 1 AU for HD 135344B. The 
molecular gas in TW Hya is found to extend inwards to 
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TABLE 2 

Best fit model parameters 



Star 


Kin [AU] 


Rout 


[AU] PA 


i 


Ah [Mq] 


SR 21 i; = 1-0 


7.6±0.4 


> 10 


16±3° 


22±40a 


O.QitO.l'' 


SR 21 v = 2-l 


7.0±0.4 


> 10 


15±4° 


20±5°^ 


l.OitO.l'" 


HD 135344B v = 


1-0 0.3±0.3 


> 15 


56±2° 


14±4oc 




TW Hya v = 1-0 


O.llitO.07 


> 1.5 


332±10° 


4.3±1.0°'^ 





''Assuming Af« = 0.7 - 1.5 Mq. 
''Assuming i = 20° (Brown 2007). 
'^Assuming M, = 1.0 - 2.0M©. 
''Assuming M, = 0.5 - 0.9 M©. 

0.1 AU, which is weU inside the 4 AU optically thin (at 
A > 1 /im) region reported by Calvet et al. (2002). 

The inclination determined from the recent Submil- 
limeter Array (SIMA) images (O'.'S resolution) of the 
SR 21 disk (Brown 2007) can be used to break the 
inclination-mass degeneracy in the spectro-astrometry 
and determine the stellar mass. Specifically, an incli- 
nation of 20° gives a stellar mass of 1.1 ±0.1 Mq for SR 
21. While HD 135344B is viewed too face-on to get a 
good estimate of the inclination from the SMA images, a 
very broad assumption on the stellar mass of 1.0-2.0 Mq, 
given the spectral type (F4), still gives a well-determined 
inclination of 14 ± 3° . Similarly, using a broad range of 
stellar masses of 0.7±0.2 Mq for TW Hya, the inclination 
of the inner disk is constrained to 4.3 ± 1.0°. 

3.2. Azimuthal asymmetries 

Because the models presented are azimuthally symmet- 
ric and assume isotropic local line emission, they produce 
symmetric (to inversion of the velocity axis) line profiles 
and symmetric (to inversion of both the velocity and the 
position axis) position-velocity diagrams. However, sig- 
nificant departures from these symmetries are apparent 
in the data, showing that the line emission from the in- 
ner disks is not azimuthally symmetric about the contin- 
uum centroid. This is especially apparent in SR 21 CO 
(u=l-0) at PA=250 and 310°, as well as in HD 135344 
at PA=240°. The spectro-astrometric technique is thus 
well-suited for detecting disk asymmetries at a specific 
azimuth angle ip. The caveat is that it cannot directly 
distinguish between an offset in radius in the line emis- 
sion at ip and an offset in the continuum emission in the 
opposite direction, + it. We define the asymmetry di- 
rection as the azimuth angle at which the line emission 
radial offset is positive, Xasym{v) = X{v) + AX{v). In 
these disks, the rotation period of the molecular emis- 
sion region is on the order of 1-10 years. A spectro- 
astrometric monitoring program should therefore be able 
to detect changes in the azimuth of the asymmetries. 

The distribution of gas in the inner disks of SR 21 
and HD 135344B can also be compared to the dust of 
the disks outside of their inner gaps, as imaged by the 
SMA (Brown 2007), clearly showing the presence of the 
dust gaps inferred from SED analyses. Fig. 4 compares 
the spectro-astrometric molecular line image information 
with these SMA images, including the asymmetry direc- 
tions. The best fit position angles as well as the asymme- 
tries appear to match the appearance of the dust disks 
very well, in particular for SR 21. The strongest submil- 
Hmeter emission blobs for both SR 21 and HD 135344B 
are aligned perpendicular to the major axis of the disk, 
indicating that the blobs are not due to a geometric ef- 




FlG. 4. — Sketch of the CO rovibrational line emission relative to 
the submillimeter continuum emission from Brown (2007) in SR 21 
and HD 135344B. The emission in HD135344B extends from 10- 
15 AU to within 0.5 AU from the central star, whereas the emission 
from SR 21 is being constrained to a relatively narrow ring around 
7AU. This is illustrated in the figure by a filled region for HD 
135344B and a ring for SR 21. The solid line indicates the best fit 
position angle and the "dumbbell" shows which side of disk is red- 
shifted; note that the rotation is measured by gas on much smaller 
size scales than the submillimeter continuum emission ring. The 
arrows indicate the line asymmetry directions, as defined in the 
text. 

feet in an optically thin disk (which would place the blobs 
parallel to the major axis), but likely indicate real tem- 
perature and/or density differences. 

4. DISCUSSION 

We have shown that spectro-astrometry of the funda- 
mental rovibrational band of CO can directly measure 
many of the basic geometric parameters (position an- 
gle, inclination, radial distribution and departures from 
axisymmetry) of molecular gas in inner protoplanetary 
disks with much fewer ambiguities than simple spec- 
troscopy. 

For the three disks studied in this paper, the high disk 
masses (Brown 2007) as well as the presence of molecular 
gas inside the dust gap argue against photo-evaporation 
as a gap-forming mechanism. However, is it possible, us- 
ing the spectro-astrometry, to distinguish between grain 
growth and dynamical clearing? Further, is it possible 
to rule out dynamical clearing by a massive, even stellar, 
companion? 

SR 21 — The presence of molecular gas in SR 21 at 7 AU 
argues against a stellar companion as the cause of the 
observed dust gap at 18 AU. However, one possible in- 
terpretation of the truncation of gas within 7AU, is the 
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presence of a relatively massive companion at < 3.5 AU, 
depending on eccentricity (Artymowicz & Lubow 1994). 
Such a companion, however, should still allow the pres- 
ence of an inner disk at 0.45 AU that may drain rapidly if 
not replenished from larger radii, as discussed in Brown 
et al. (2007). 

Could the dust gap between 7 and 18 AU be due to 
grain growth, but with the gas continuing inwards as a 
natural extension of the outer disk? By definition, dust 
in protoplanetary disk gaps and holes has much lower 
optical depth in the infrared in the vertical direction 
than disks without such structure. In this case, infrared 
rovibrational lines will be much better tracers of the to- 
tal vertical gas column density than normally, since the 
disk midplanes will no longer be blocked from view by 
a high optical depth due to dust continuum absorption. 
Indeed, the presence of strong lines of rare isotopologues 
(i^CO, C^^O and C"0, sec Figure 1) in emission in SR 
21 indicates higher vertical gas column densities than HD 
135344B; SR 21 has ^^CO/^CO line ratios of ~ 3, while 
HD 135344B has ratios of ~ 8 - 10. A high vertical gas 
column density in SR 21 argues for dust gap formation 
by grain growth beyond 7 AU. A truncation of molecular 
gas within 7AU, however, is more consistent with dy- 
namical clearing by a massive companion at ~ 3 — 4 AU. 
It thus seems that more than one mechanism is in play in 
the SR 21 disk. More detailed multi-dimensional radia- 
tive transfer modeling is needed to accurately determine 
the vertical CO column densities. 

HD I35344B — Our detection of molecular gas extending 
from within 0.3 AU of the star, to at least 15 AU argues 
against clearing by a sfe//ar companion for HD 135344B, 
since this would likely have led to a completely evacu- 
ated inner hole. Further, the vertical column density of 
molecular gas seems relatively low, as evidenced by the 
low i2co/"CO line ratios. Since HD 135344B is sur- 
rounded by an otherwise massive disk, the CO rovibra- 
tional lines are more consistent with it having an actual 
density maximum at the outer dust gap edge, possibly 
caused by a planetary mass companion orbiting at 10- 
20 AU. In this case, care must be taken to rule out that 
photo-dissociation has lowered the vertical column den- 
sity of CO, which might happen if grain growth has re- 
moved the UV dust opacity. If this happens, molecular 
gas will no longer be a good tracer of the total gas col- 
umn. However, even in the absence of opacity-producing 
dust, CO will generally become self shielding, except for 
very low mass disks (Jonkheid et al. 2006). 

TW Hya — Our observed distribution of warm CO gas 
between 0.11 - 1.5 AU for TW Hya shows that molecular 
gas is present in the optically thin dust region found by 
interferometric dust measurements. It is also consistent 
with the recent detection of a young planet at 0.04 AU, 
which is likely to truncate the disk at ~0.1AU (Seti- 
awan et al. 2008). The i^co/^^co line ratio in TW 
Hya is very high (> 15, see also Salyk et al. 2007), in- 
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dicating a vertical gas column density that is lower than 
expected from a smooth, monotonic radial surface den- 
sity profile. Thus the data indicate that a mechanism 
other than grain growth is producing the optically thin 
inner disk hole. Note that if the 4 AU inner hole is being 
dynamically cleared by a planetary companion, it will 
not be the companion detected at 0.04 AU, but another 
one with an orbit at 1-2 AU. 

Our estimate of the inclination of the inner disk for 
TW Hya of 4.3 ± 1.0° brings the TW Hya companion 
mass to 16^3 Mjup. The inner disk inclination, if more 
representative of the companion orbit inclination than 
the outer disk inclination of 7 ± 1° by Qi et al. (2004), 
places this short period companion in the "brown dwarf 
desert" (Grether & Lineweaver 2006). Further, the dif- 
ference in inclination measured in the inner and outer 
disks of TW Hya is indicative of a warp in the disk, also 
consistent with the presence of a perturbing young plan- 
etary system. A warp in a protoplanetary disk can be 
induced by an encounter with a passing star, or through 
dynamical interactions with a planetary system such as 
has been suggested for the famous Beta Pictoris debris 
disk (Golimowski et al. 2006). There are no obvious 
nearby cluster members to the TW Hya disk, so a stel- 
lar encounter is a less likely explanation for a disk warp. 
A disk warp induced by a planet probably requires that 
the inclination of the orbit of the perturbing planet has 
been changed relative to the mean disk inclination by an 
interaction with another planet. It is therefore tempting 
to speculate that there is an additional planet in the TW 
Hya system, possibly causing the observed 4AU gap. 

This paper presents the first results of an extensive 
spectro-astrometric imaging survey of molecular gas in 
the planet forming regions of nearby protoplanetary 
disks. With the powerful resolution of CRIRES of < 1 
milli-arcsec, comparable to or even better than the reso- 
lution of near-infrared interferometers, the CO rovibra- 
tional emission from nearby protoplanetary disks is gen- 
erally resolved both spatially and kinematically. A sen- 
sitive interferometric capability at 3-5 /xm with resolv- 
ing powers of 50,000-100,000 will be required to improve 
on the spectro-astromctry presented here. Spatially re- 
solving the emission from other molecules, in particular 
the water lines recently reported in Salyk et al. (2008), 
is likely within the capability of CRIRES, and would 
directly probe the chemical variations in the planet- 
forming region of circumstellar disks. 
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